Alternative epidemic indicators for COVID-19 in three settings with incomplete death registration systems.

Not all COVID-19 deaths are officially reported, and particularly in low-income and humanitarian settings, the magnitude of reporting gaps remains sparsely characterized. Alternative data sources, including burial site worker reports, satellite imagery of cemeteries, and social media-conducted surveys of infection may offer solutions. By merging these data with independently conducted, representative serological studies within a mathematical modeling framework, we aim to better understand the range of underreporting using examples from three major cities: Addis Ababa (Ethiopia), Aden (Yemen), and Khartoum (Sudan) during 2020. We estimate that 69 to 100%, 0.8 to 8.0%, and 3.0 to 6.0% of COVID-19 deaths were reported in each setting, respectively. In future epidemics, and in settings where vital registration systems are limited, using multiple alternative data sources could provide critically needed, improved estimates of epidemic impact. However, ultimately, these systems are needed to ensure that, in contrast to COVID-19, the impact of future pandemics or other drivers of mortality is reported and understood worldwide.

Lessons from COVID-19 for rescalable data collection.

Novel data and analyses have had an important role in informing the public health response to the COVID-19 pandemic. Existing surveillance systems were scaled up, and in some instances new systems were developed to meet the challenges posed by the magnitude of the pandemic. We describe the routine and novel data that were used to address urgent public health questions during the pandemic, underscore the challenges in sustainability and equity in data generation, and highlight key lessons learnt for designing scalable data collection systems to support decision making during a public health crisis. As countries emerge from the acute phase of the pandemic, COVID-19 surveillance systems are being scaled down. However, SARS-CoV-2 resurgence remains a threat to global health security; therefore, a minimal cost-effective system needs to remain active that can be rapidly scaled up if necessary. We propose that a retrospective evaluation to identify the cost-benefit profile of the various data streams collected during the pandemic should be on the scientific research agenda.

Understanding the Potential Impact of Different Drug Properties on Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Transmission and Disease Burden: A Modelling Analysis.

BACKGROUND: The public health impact of the coronavirus disease 2019 (COVID-19) pandemic has motivated a rapid search for potential therapeutics, with some key successes. However, the potential impact of different treatments, and consequently research and procurement priorities, have not been clear. METHODS: Using a mathematical model of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) transmission, COVID-19 disease and clinical care, we explore the public-health impact of different potential therapeutics, under a range of scenarios varying healthcare capacity, epidemic trajectories; and drug efficacy in the absence of supportive care. RESULTS: The impact of drugs like dexamethasone (delivered to the most critically-ill in hospital and whose therapeutic benefit is expected to depend on the availability of supportive care such as oxygen and mechanical ventilation) is likely to be limited in settings where healthcare capacity is lowest or where uncontrolled epidemics result in hospitals being overwhelmed. As such, it may avert 22% of deaths in high-income countries but only 8% in low-income countries (assuming R = 1.35). Therapeutics for different patient populations (those not in hospital, early in the course of infection) and types of benefit (reducing disease severity or infectiousness, preventing hospitalization) could have much greater benefits, particularly in resource-poor settings facing large epidemics. CONCLUSIONS: Advances in the treatment of COVID-19 to date have been focused on hospitalized-patients and predicated on an assumption of adequate access to supportive care. Therapeutics delivered earlier in the course of infection that reduce the need for healthcare or reduce infectiousness could have significant impact, and research into their efficacy and means of delivery should be a priority.

Estimating the COVID-19 infection fatality ratio accounting for seroreversion using statistical modelling.

Background: The infection fatality ratio (IFR) is a key statistic for estimating the burden of coronavirus disease 2019 (COVID-19) and has been continuously debated throughout the COVID-19 pandemic. The age-specific IFR can be quantified using antibody surveys to estimate total infections, but requires consideration of delay-distributions from time from infection to seroconversion, time to death, and time to seroreversion (i.e. antibody waning) alongside serologic test sensitivity and specificity. Previous IFR estimates have not fully propagated uncertainty or accounted for these potential biases, particularly seroreversion. Methods: We built a Bayesian statistical model that incorporates these factors and applied this model to simulated data and 10 serologic studies from different countries. Results: We demonstrate that seroreversion becomes a crucial factor as time accrues but is less important during first-wave, short-term dynamics. We additionally show that disaggregating surveys by regions with higher versus lower disease burden can inform serologic test specificity estimates. The overall IFR in each setting was estimated at 0.49-2.53%. Conclusion: We developed a robust statistical framework to account for full uncertainties in the parameters determining IFR. We provide code for others to apply these methods to further datasets and future epidemics.

A framework for reconstructing SARS-CoV-2 transmission dynamics using excess mortality data.

The transmission dynamics and burden of SARS-CoV-2 in many regions of the world is still largely unknown due to the scarcity of epidemiological analyses and lack of testing to assess the prevalence of disease. In this work, we develop a quantitative framework based on excess mortality data to reconstruct SARS-CoV-2 transmission dynamics and assess the level of underreporting in infections and deaths. Using weekly all-cause mortality data from Iran, we are able to show a strong agreement between our attack rate estimates and seroprevalence measurements in each province and find significant heterogeneity in the level of exposure across the country with 11 provinces reaching near 100% attack rates. Despite having a young population, our analysis reveals that incorporating limited access to medical services in our model, coupled with undercounting of COVID-19-related deaths, leads to estimates of infection fatality rate in most provinces of Iran that are comparable to high-income countries.

Modelling the impact of vaccine hesitancy in prolonging the need for Non-Pharmaceutical Interventions to control the COVID-19 pandemic.

Background: Vaccine hesitancy - a delay in acceptance or refusal of vaccines despite availability - has the potential to threaten the successful roll-out of SARS-CoV-2 vaccines globally. In this study, we aim to understand the likely impact of vaccine hesitancy on the control of the COVID-19 pandemic. Methods: We modelled the potential impact of vaccine hesitancy on the control of the pandemic and the relaxation of non-pharmaceutical interventions (NPIs) by combining an epidemiological model of SARS-CoV-2 transmission with data on vaccine hesitancy from population surveys. Results: Our simulations suggest that the mortality over a 2-year period could be up to 7.6 times higher in countries with high vaccine hesitancy compared to an ideal vaccination uptake if NPIs are relaxed. Alternatively, high vaccine hesitancy could prolong the need for NPIs to remain in place. Conclusions: While vaccination is an individual choice, vaccine-hesitant individuals have a substantial impact on the pandemic trajectory, which may challenge current efforts to control COVID-19. In order to prevent such outcomes, addressing vaccine hesitancy with behavioural interventions is an important priority in the control of the COVID-19 pandemic.

Reassessing Reported Deaths and Estimated Infection Attack Rate during the First 6 Months of the COVID-19 Epidemic, Delhi, India.

India reported >10 million coronavirus disease (COVID-19) cases and 149,000 deaths in 2020. To reassess reported deaths and estimate incidence rates during the first 6 months of the epidemic, we used a severe acute respiratory syndrome coronavirus 2 transmission model fit to data from 3 serosurveys in Delhi and time-series documentation of reported deaths. We estimated 48.7% (95% credible interval 22.1%-76.8%) cumulative infection in the population through the end of September 2020. Using an age-adjusted overall infection fatality ratio based on age-specific estimates from mostly high-income countries, we estimated that just 15.0% (95% credible interval 9.3%-34.0%) of COVID-19 deaths had been reported, indicating either substantial underreporting or lower age-specific infection-fatality ratios in India than in high-income countries. Despite the estimated high attack rate, additional epidemic waves occurred in late 2020 and April-May 2021. Future dynamics will depend on the duration of natural and vaccine-induced immunity and their effectiveness against new variants.

Communicating uncertainty in epidemic models.

While mathematical models of disease transmission are widely used to inform public health decision-makers globally, the uncertainty inherent in results are often poorly communicated. We outline some potential sources of uncertainty in epidemic models, present traditional methods used to illustrate uncertainty and discuss alternative presentation formats used by modelling groups throughout the COVID-19 pandemic. Then, by drawing on the experience of our own recent modelling, we seek to contribute to the ongoing discussion of how to improve upon traditional methods used to visualise uncertainty by providing a suggestion of how this can be presented in a clear and simple manner.

Reproducible parallel inference and simulation of stochastic state space models using odin, dust, and mcstate.

State space models, including compartmental models, are used to model physical, biological and social phenomena in a broad range of scientific fields. A common way of representing the underlying processes in these models is as a system of stochastic processes which can be simulated forwards in time. Inference of model parameters based on observed time-series data can then be performed using sequential Monte Carlo techniques. However, using these methods for routine inference problems can be made difficult due to various engineering considerations: allowing model design to change in response to new data and ideas, writing model code which is highly performant, and incorporating all of this with up-to-date statistical techniques. Here, we describe a suite of packages in the R programming language designed to streamline the design and deployment of state space models, targeted at infectious disease modellers but suitable for other domains. Users describe their model in a familiar domain-specific language, which is converted into parallelised C++ code. A fast, parallel, reproducible random number generator is then used to run large numbers of model simulations in an efficient manner. We also provide standard inference and prediction routines, though the model simulator can be used directly if these do not meet the user's needs. These packages provide guarantees on reproducibility and performance, allowing the user to focus on the model itself, rather than the underlying computation. The ability to automatically generate high-performance code that would be tedious and time-consuming to write and verify manually, particularly when adding further structure to compartments, is crucial for infectious disease modellers. Our packages have been critical to the development cycle of our ongoing real-time modelling efforts in the COVID-19 pandemic, and have the potential to do the same for models used in a number of different domains.

Country differences in transmissibility, age distribution and case-fatality of SARS-CoV-2: a global ecological analysis.

Objectives The first COVID-19 pandemic waves in many low-income countries appeared milder than initially forecasted. We conducted a country-level ecological study to describe patterns in key SARS-CoV-2 outcomes by country and region and explore associations with potential explanatory factors, including population age structure and prior exposure to endemic parasitic infections. Methods We collected publicly available data and compared them using standardisation techniques. We then explored the association between exposures and outcomes using random forest and linear regression. We adjusted for potential confounders and plausible effect modifications. Results While mean time-varying reproduction number was highest in the European and Americas regions, median age of death was lower in the Africa region, with a broadly similar case-fatality ratio. Population age was strongly associated with mean (ß=0.01, 95% CI, 0.005, 0.011) and median age of cases (ß=-0.40, 95% CI, -0.53, -0.26) and deaths (ß= 0.40, 95% CI, 0.17, 0.62). Conclusions Population age seems an important country-level factor explaining both transmissibility and age distribution of observed cases and deaths. Endemic infections seem unlikely, from this analysis, to be key drivers of the variation in observed epidemic trends. Our study was limited by the availability of outcome data and its causally uncertain ecological design.

Key epidemiological drivers and impact of interventions in the 2020 SARS-CoV-2 epidemic in England.

We fitted a model of SARS-CoV-2 transmission in care homes and the community to regional surveillance data for England. Compared with other approaches, our model provides a synthesis of multiple surveillance data streams into a single coherent modeling framework, allowing transmission and severity to be disentangled from features of the surveillance system. Of the control measures implemented, only national lockdown brought the reproduction number (Rt eff) below 1 consistently; if introduced 1 week earlier, it could have reduced deaths in the first wave from an estimated 48,600 to 25,600 [95% credible interval (CrI): 15,900 to 38,400]. The infection fatality ratio decreased from 1.00% (95% CrI: 0.85 to 1.21%) to 0.79% (95% CrI: 0.63 to 0.99%), suggesting improved clinical care. The infection fatality ratio was higher in the elderly residing in care homes (23.3%, 95% CrI: 14.7 to 35.2%) than those residing in the community (7.9%, 95% CrI: 5.9 to 10.3%). On 2 December 2020, England was still far from herd immunity, with regional cumulative infection incidence between 7.6% (95% CrI: 5.4 to 10.2%) and 22.3% (95% CrI: 19.4 to 25.4%) of the population. Therefore, any vaccination campaign will need to achieve high coverage and a high degree of protection in vaccinated individuals to allow nonpharmaceutical interventions to be lifted without a resurgence of transmission.

Using syndromic measures of mortality to capture the dynamics of COVID-19 in Java, Indonesia, in the context of vaccination rollout.

BACKGROUND: As in many countries, quantifying COVID-19 spread in Indonesia remains challenging due to testing limitations. In Java, non-pharmaceutical interventions (NPIs) were implemented throughout 2020. However, as a vaccination campaign launches, cases and deaths are rising across the island. METHODS: We used modelling to explore the extent to which data on burials in Jakarta using strict COVID-19 protocols (C19P) provide additional insight into the transmissibility of the disease, epidemic trajectory, and the impact of NPIs. We assess how implementation of NPIs in early 2021 will shape the epidemic during the period of likely vaccine rollout. RESULTS: C19P burial data in Jakarta suggest a death toll approximately 3.3 times higher than reported. Transmission estimates using these data suggest earlier, larger, and more sustained impact of NPIs. Measures to reduce sub-national spread, particularly during Ramadan, substantially mitigated spread to more vulnerable rural areas. Given current trajectory, daily cases and deaths are likely to increase in most regions as the vaccine is rolled out. Transmission may peak in early 2021 in Jakarta if current levels of control are maintained. However, relaxation of control measures is likely to lead to a subsequent resurgence in the absence of an effective vaccination campaign. CONCLUSIONS: Syndromic measures of mortality provide a more complete picture of COVID-19 severity upon which to base decision-making. The high potential impact of the vaccine in Java is attributable to reductions in transmission to date and dependent on these being maintained. Increases in control in the relatively short-term will likely yield large, synergistic increases in vaccine impact.

Temperature and population density influence SARS-CoV-2 transmission in the absence of nonpharmaceutical interventions.

As COVID-19 continues to spread across the world, it is increasingly important to understand the factors that influence its transmission. Seasonal variation driven by responses to changing environment has been shown to affect the transmission intensity of several coronaviruses. However, the impact of the environment on severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) remains largely unknown, and thus seasonal variation remains a source of uncertainty in forecasts of SARS-CoV-2 transmission. Here we address this issue by assessing the association of temperature, humidity, ultraviolet radiation, and population density with estimates of transmission rate (R). Using data from the United States, we explore correlates of transmission across US states using comparative regression and integrative epidemiological modeling. We find that policy intervention ("lockdown") and reductions in individuals' mobility are the major predictors of SARS-CoV-2 transmission rates, but, in their absence, lower temperatures and higher population densities are correlated with increased SARS-CoV-2 transmission. Our results show that summer weather cannot be considered a substitute for mitigation policies, but that lower autumn and winter temperatures may lead to an increase in transmission intensity in the absence of policy interventions or behavioral changes. We outline how this information may improve the forecasting of COVID-19, reveal its future seasonal dynamics, and inform intervention policies.

Leveraging community mortality indicators to infer COVID-19 mortality and transmission dynamics in Damascus, Syria.

The COVID-19 pandemic has resulted in substantial mortality worldwide. However, to date, countries in the Middle East and Africa have reported considerably lower mortality rates than in Europe and the Americas. Motivated by reports of an overwhelmed health system, we estimate the likely under-ascertainment of COVID-19 mortality in Damascus, Syria. Using all-cause mortality data, we fit a mathematical model of COVID-19 transmission to reported mortality, estimating that 1.25% of COVID-19 deaths (sensitivity range 1.00% - 3.00%) have been reported as of 2 September 2020. By 2 September, we estimate that 4,380 (95% CI: 3,250 - 5,550) COVID-19 deaths in Damascus may have been missed, with 39.0% (95% CI: 32.5% - 45.0%) of the population in Damascus estimated to have been infected. Accounting for under-ascertainment corroborates reports of exceeded hospital bed capacity and is validated by community-uploaded obituary notifications, which confirm extensive unreported mortality in Damascus.

Genetic evidence for the association between COVID-19 epidemic severity and timing of non-pharmaceutical interventions.

Unprecedented public health interventions including travel restrictions and national lockdowns have been implemented to stem the COVID-19 epidemic, but the effectiveness of non-pharmaceutical interventions is still debated. We carried out a phylogenetic analysis of more than 29,000 publicly available whole genome SARS-CoV-2 sequences from 57 locations to estimate the time that the epidemic originated in different places. These estimates were examined in relation to the dates of the most stringent interventions in each location as well as to the number of cumulative COVID-19 deaths and phylodynamic estimates of epidemic size. Here we report that the time elapsed between epidemic origin and maximum intervention is associated with different measures of epidemic severity and explains 11% of the variance in reported deaths one month after the most stringent intervention. Locations where strong non-pharmaceutical interventions were implemented earlier experienced much less severe COVID-19 morbidity and mortality during the period of study.

Within-country age-based prioritisation, global allocation, and public health impact of a vaccine against SARS-CoV-2: A mathematical modelling analysis.

The worldwide endeavour to develop safe and effective COVID-19 vaccines has been extraordinary, and vaccination is now underway in many countries. However, the doses available in 2021 are likely to be limited. We extend a mathematical model of SARS-CoV-2 transmission across different country settings to evaluate the public health impact of potential vaccines using WHO-developed target product profiles. We identify optimal vaccine allocation strategies within- and between-countries to maximise averted deaths under constraints on dose supply. We find that the health impact of SARS-CoV-2 vaccination depends on the cumulative population-level infection incidence when vaccination begins, the duration of natural immunity, the trajectory of the epidemic prior to vaccination, and the level of healthcare available to effectively treat those with disease. Within a country we find that for a limited supply (doses for < 20% of the population) the optimal strategy is to target the elderly. However, with a larger supply, if vaccination can occur while other interventions are maintained, the optimal strategy switches to targeting key transmitters to indirectly protect the vulnerable. As supply increases, vaccines that reduce or block infection have a greater impact than those that prevent disease alone due to the indirect protection provided to high-risk groups. Given a 2 billion global dose supply in 2021, we find that a strategy in which doses are allocated to countries proportional to population size is close to optimal in averting deaths and aligns with the ethical principles agreed in pandemic preparedness planning.

Modelling intensive care unit capacity under different epidemiological scenarios of the COVID-19 pandemic in three Western European countries.

BACKGROUND: The coronavirus disease 2019 (COVID-19) pandemic has placed enormous strain on intensive care units (ICUs) in Europe. Ensuring access to care, irrespective of COVID-19 status, in winter 2020-2021 is essential. METHODS: An integrated model of hospital capacity planning and epidemiological projections of COVID-19 patients is used to estimate the demand for and resultant spare capacity of ICU beds, staff and ventilators under different epidemic scenarios in France, Germany and Italy across the 2020-2021 winter period. The effect of implementing lockdowns triggered by different numbers of COVID-19 patients in ICUs under varying levels of effectiveness is examined, using a 'dual-demand' (COVID-19 and non-COVID-19) patient model. RESULTS: Without sufficient mitigation, we estimate that COVID-19 ICU patient numbers will exceed those seen in the first peak, resulting in substantial capacity deficits, with beds being consistently found to be the most constrained resource. Reactive lockdowns could lead to large improvements in ICU capacity during the winter season, with pressure being most effectively alleviated when lockdown is triggered early and sustained under a higher level of suppression. The success of such interventions also depends on baseline bed numbers and average non-COVID-19 patient occupancy. CONCLUSION: Reductions in capacity deficits under different scenarios must be weighed against the feasibility and drawbacks of further lockdowns. Careful, continuous decision-making by national policymakers will be required across the winter period 2020-2021.

Reduction in mobility and COVID-19 transmission.

In response to the COVID-19 pandemic, countries have sought to control SARS-CoV-2 transmission by restricting population movement through social distancing interventions, thus reducing the number of contacts. Mobility data represent an important proxy measure of social distancing, and here, we characterise the relationship between transmission and mobility for 52 countries around the world. Transmission significantly decreased with the initial reduction in mobility in 73% of the countries analysed, but we found evidence of decoupling of transmission and mobility following the relaxation of strict control measures for 80% of countries. For the majority of countries, mobility explained a substantial proportion of the variation in transmissibility (median adjusted R-squared: 48%, interquartile range - IQR - across countries [27-77%]). Where a change in the relationship occurred, predictive ability decreased after the relaxation; from a median adjusted R-squared of 74% (IQR across countries [49-91%]) pre-relaxation, to a median adjusted R-squared of 30% (IQR across countries [12-48%]) post-relaxation. In countries with a clear relationship between mobility and transmission both before and after strict control measures were relaxed, mobility was associated with lower transmission rates after control measures were relaxed indicating that the beneficial effects of ongoing social distancing behaviours were substantial.

Evidence of initial success for China exiting COVID-19 social distancing policy after achieving containment.

Background: The COVID-19 epidemic was declared a Global Pandemic by WHO on 11 March 2020. By 24 March 2020, over 440,000 cases and almost 20,000 deaths had been reported worldwide. In response to the fast-growing epidemic, which began in the Chinese city of Wuhan, Hubei, China imposed strict social distancing in Wuhan on 23 January 2020 followed closely by similar measures in other provinces. These interventions have impacted economic productivity in China, and the ability of the Chinese economy to resume without restarting the epidemic was not clear. Methods: Using daily reported cases from mainland China and Hong Kong SAR, we estimated transmissibility over time and compared it to daily within-city movement, as a proxy for economic activity. Results: Initially, within-city movement and transmission were very strongly correlated in the five mainland provinces most affected by the epidemic and Beijing. However, that correlation decreased rapidly after the initial sharp fall in transmissibility. In general, towards the end of the study period, the correlation was no longer apparent, despite substantial increases in within-city movement. A similar analysis for Hong Kong shows that intermediate levels of local activity were maintained while avoiding a large outbreak. At the very end of the study period, when China began to experience the re-introduction of a small number of cases from Europe and the United States, there is an apparent up-tick in transmission. Conclusions: Although these results do not preclude future substantial increases in incidence, they suggest that after very intense social distancing (which resulted in containment), China successfully exited its lockdown to some degree. Elsewhere, movement data are being used as proxies for economic activity to assess the impact of interventions. The results presented here illustrate how the eventual decorrelation between transmission and movement is likely a key feature of successful COVID-19 exit strategies.